security  ci*>u.icatiom  or  tnisraoe  m>iw  omm  lnf«M 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


IRiEnT'S  CATALOG  NIMH* 


«  TITUS  <m*  luMMtA) 

A  FURTHER  COMPARISON  OF  DETERMINISTIC  AND 
STOCHASTIC  LANCHESTER-TYPE  COMBAT  MODELS 


T.  Au  ThO HCI 


i.  tyre  or  mbromt  a  remoo  covered 

Master's  Thesis; 

September  1980 


ft.  PinroitMiMa  ona.  report  number 


Kurt  Dieter  Klemm 


RpfpMRf  OmRURTAM  NAME  AMO  AOOftCIft 

Naval  Postgraduate  School 
Monterey,  California  93940 


II  CONTROLLING  ORRlCE  WANS  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


rmooran  element.  rmoject.  task 
ANSA  A  ROMS  UNIT  HUHOSNS 


IS.  REPORT  OATS 

September  1980 


f«.  MON>TO«»NG  AGCNCV  nam«  *  lOQtfllflf  tff/ferwr t  frmm  CmmtrmiHmj  OtH *•)  1ft.  |ICU*ITV  CLASS,  (mi  t*tm  MRaM) 

Naval  Postgraduate  School  Unclassified 

Monterey,  California  93940  ,,  .  _ _ 

IE  Jj^C^ttiriCATION/OORNONAOlHC 


IS.  OlSTMl SuTlON  STATSNSMT  Cl  MU*  Mnra.ii 

Approved  for  public  release;  distribution  unlimited 


IT.  DISTRIBUTION  ST  ATSNSNT  Cl  »•  AA.'.AAI  «NnN  M  IlNl  M,  II  RNmM  I 


IS  KCY  ROMOS  1C ahiRma  AH  HWH  XAa  1/  a ac.n .ait  anA  iRsnrllr  Sr  Sirs*  rrAr) 


Lanchester,  Square-law  Attrition,  Stochastic,  Combat  Models,  Bivariate  Death 
Process 


SON*  ABSTRACT  rCMlfUM  MR  mntn  «<#•  II  ammamt  MR#  IBmi|R>  If  Urn 

This  thesis  examines  the  differences  of  deterministic  and  stochastic 
LANCHESTER-type  combat  models.  Using  an  example  of  square-law  attrition, 
solution  methods  and  solutions  are  described.  A  new  analytic  solution  for 
equal  attrition  rate  coefficients  is  given.  The  numerical  comparison  includes 
hypotheses  about  the  expected  force  levels  and  the  variability  in  the  expected 
force  levels  as  a  function  of  time,  initial  force  levels,  and  breakpoint  force 
levels. 


00  ,sr»  1<73  EDITION  OR  '  NOV  ••  It  OBSOLETE 

(Page  1)  S/N  0  10  I- 01  4- MO  I 


IECUMITT  CL AftttRICATION  OR  THIE  MADE  I 


Approved  for  public  release;  distribution  unlimited 


A  Further  Comparison  of  Deterministic  and  Stochastic 
Lanchester-Type  Combat  Models 


by 

Kurt  Dieter  Klemm 

Captain,  Army,  Federal  Republic  of  Germany 
Ingenieur  (grad.)»  Fachhochschule  des  Heeres  1,  1975 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 
from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1980 


Author 
Approved  by: 


Second  Reader 


Chairman,  Department  bf  Operatipns  Research 

A/ 

Dean  of  Information  and  Policy  Sciences 


Z 


ABSTRACT 


This  thesis  examines  the  differences  of  deterministic  and  stochastic 

LANCHESTER-type  combat  models.  Using  an  example  of  square-law  attrition, 
solution  methods  and  solutions  are  described.  A  new  analytic  solution 
for  equal  attrition  rate  coefficients  is  given.  The  numerical  comparison 
includes  hypotheses  about  the  expected  force  levels  and  the  variability 

in  the  expected  force  levels  as  a  function  of  time,  initial  force  levels, 

and  breakpoint  force  levels. 
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I.  INTRODUCTION 


Combat  models  are  widely  used  as  decision  aids  in  the  defense¬ 
planning  process,  at  least  within  the  NATO  alliance.  Current  operation¬ 
al  combat  models  are  very  complex  because  combat  is  a  very  complex 
process.  Unfortunately  it  is  difficult  (if  not  impossible)  for  the 
beginner  to  understand  the  modelling  approaches,  concepts  and  motivation, 
that  may  have  been  used  to  build  such  operational  models.  However,  one 
frequently  considers  a  simple  model  as  a  paradigm  for  the  development  and 
understanding  of  such  complex  models.  This  basic  approach  will  be  used 
in  this  thesis  to  explore  certain  issues  in  the  on-going  debate  about 
the  relative  merits  of  stochastic  and  deterministic  combat  models. 

A  simple  model  is  examined  to  explore  differences  between  a  determin¬ 
istic  and  a  stochastic  approach  to  a  certain  type  of  analytical  combat 
model.  As  already  mentioned,  combat  is  a  very  complex  process,  but  it 
is  also  a  complex  random  process,  which  can  be  supported  by  many  examples 
from  military  history.  Analytical  models  are  abstractions  and  very  often 
simplifications  of  reality.  It  seems  to  be  a  legitimate  question  to  ask, 
what  effects  the  further  abstraction  of  neglecting  the  randomness  in 
combat  may  have.  At  this  moment,  it  should  be  pointed  out  that  within 
existing  operational  analytical  models,  both  stochastic  and  deterministic 
models  are  used. 

Previous  work  done  by  SPRINGALL  [9]  and  CLARK  [4]  evolved  around 
theoretical  aspects.  Their  main  concern  was  to  give  exact  analytical 
solutions  and  their  proofs.  CRAIG  [5]  started  to  explore  the  differences 
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between  stochastic  and  deterministic  models  more  from  the  numerical 
point  of  view,  which  will  be  continued  in  this  thesis. 

In  the  next  chapter,  a  deterministic  and  stochastic  version  of  a 
differential  combat  model  will  be  described.  The  deterministic  versions 
are  well-known  as  LANCHESTER's  equations  of  modern  warfare,  which  were 
developed  in  1914.  Combat  models,  which  model  attrition  from  enemy 
action  through  a  system  of  differential  equations,  are  usually  referred 
to  as  LANCHESTER-type  models  of  warfare. 
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I I .  THE  PARADIGM  MODELS 


A.  THE  DETERMINISTIC  MODEL 

First,  LANCHESTER's  equations  of  modern  warfare  will  be  briefly 
reviewed  and  some  simple  extensions  given. 

In  1914  LANCHESTER  [7]  hypothesized  that  under  "modern  conditions" 
in  a  combat  between  two  homogeneous  forces  the  firepower  of  the  surviving 
weapons  of  one  side  can  be  concentrated  on  the  surviving  targets  of  the 
other  side,  so  that  each  side's  casualty  rate  is  proportional  to  the  num¬ 
ber  of  enemy  firers.  This  can  be  described  by  the  following  equations: 


at*  -w 

(2.1) 

(2.2) 

i  initial  conditions 

O 

X 

II 

o 

X 

(2.3) 

o 

>> 

II 

o 

>- 

(2.4) 

where  a  is  the  attrition  rate  with  which  the  Y-force  attrits  the 
X-force,  similarly  for  b.  Xg  and  Yg  are  the  initial  force  levels  and 
X ( t )  and  Y(t)  are  the  force  levels  at  time  t.  The  force  levels,  as  a 
function  of  time  t,  can  be  written  as 


X(t)=XgCOsh(Vab  t)-Va7b  ygSinhO/ab’ t) 

(2.5) 

Y(t)=ygCOsh(<yab’  t)-Vb7a  XgSinh(Vab  t) 

(2.6) 

The  state  equation  relating  initial  force  levels  with  force  levels  at 
some  time  t  can  be  derived  by  dividing  (2.1)  by  (2.2),  separating 
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variables  and  integrating  to  yield 

b(x02-X(t)2)=a(y02-Y(t)2).  (2.7) 

This  form  of  the  state  equation  explains  why  this  model  is  also  referred 
to  as  the  "square-law"  attrition  process.  WEISS  [11]  has  given  a  set  of 
assumptions  under  which  LANCHESTER's  equations  for  modern  warfare  may 
apply: 

Al.)  Two  homogeneous  forces  are  engaged  in  combat.  Every  unit 
on  a  particular  side  has  the  same  capabilities.  The  attrition 
rate  may  be  different  for  the  two  forces. 

A2.)  Each  unit  on  one  side  is  within  weapon  range  of  all  units 
on  the  other  side. 

A3.)  The  effects  of  successive  rounds  on  the  target  are 
independent. 

A4.)  Each  unit  has  perfect  knowledge  of  target  locations  and 
fires  only  at  live  target  (one  at  a  time)  killing  them  at  a 
constant  rate,  which  does  not  depend  on  the  number  of  targets 
alive. 

A5.)  Fire  is  uniformly  distributed  over  surviving  targets. 

The  above  model  implies  a  fight  until  one  force  is  annihilated. 
Therefore  the  model  will  be  slightly  changed  by  introducing  the  concept 
of  unit  breakpoints,  Xfap  and  Ybp,  which  are  simply  force  levels  at  which 
the  side,  who  reaches  it  first  "breaks  off"  the  engagement,  leaving  the 
other  side  as  a  winner.  Also,  to  be  more  precise,  it  should  be  noted 
that  negative  force  levels  for  breakpoints  equal  zero  or  force  levels 
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less  than  nonzero  breakpoints  are  impossible.  So  the  deterministic 
LANCHESTER-type  combat  model  with  "square-law"  attrition  takes  the 
following  form: 

-ay  xbp<X(t) 

ybp<Y(t)  (2.8) 

0  otherwise 


£= 


-bx  xbp<X(t) 

ybp<Y<t> 

0  otherwise 


with  initial  conditions 
X(0)  =  xQ 
Y(0)  -  yQ 


(2.9) 


(2.10) 

(2.11) 


The  model  in  this  form,  equations  (2.8)  through  (2.11),  was  used  for 
comparisons  throughout  the  thesis. 


B.  THE  STOCHASTIC  MODEL 

There  are  several  ways  to  include  random  variations  in  LANCHESTER- 
type  models.  These  are: 

*  The  attrition  rate  coefficients  may  be  random  variables. 

*  The  enemy's  initial  force  level  may  be  a  random  variable, 
weakening  the  assumption  of  perfect  knowledge. 

*  The  breakpoints  may  be  random  variables. 

*  The  casualty  rate  is  fixed,  but  the  occurrence  of  casualties 
over  time  may  be  random. 

The  only  random  variation  considered  here  will  be  the  random  occur¬ 
rence  of  casualties  over  time.  Another  specification  was  to  choose  a 
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model  si> '"Jar  to  the  "square-1 aw"attriti on  in  order  to  allow  comparisons. 

In  other  words,  the  question  to  be  asked  is  "how  do  random  fluctuations 
in  the  occurrence  of  casualties  modify  the  deterministic  results  of  the 
square-law  attrition  process?" 

The  approach  used  here  was  a  continuous  parameter  MARKOV  chain  model, 
where  the  time  t  varies  continuously  and  the  number  of  combatants  on 
each  side  is  a  non-negative  integer.  Let  M(t)  be  the  size  of  the  X-force 
at  time  t  with  a  particular  state  value  m.  Let  N(t)  be  the  size  of  the 
Y-force  at  time  t  with  a  particular  state  value  n.  Let  mQ  and  nQ  be  the 
initial  force  levels  and  m^  ,  n^  be  the  breakpoint  force  levels  of  X 
and  Y  respectively.  Fig.  1  shows  the  state  space  of  this  MARKOV  chain 
model.  Note  that  at  a  given  time  t,  any  state  is  described  by  the  two 
force  levels  of  the  X  and  Y  force.  As  each  side  loses  units  due  to 
attrition  and  no  replacements  are  allowed,  it  is  easy  to  understand  why 
BILLARD  [1]  referred  to  this  type  of  process  as  a  "bivariate  death  process." 
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For  the  description  of  the  stochastic  square-law  attrition  process 
corresponding  to  the  two  deterministic  differential  equations  (2.8)  and 
(2.9),  a  system  of  many  differential  equations,  depending  on  the  battle 
termination  model,  is  required.  This  system  will  be  given  for  a  fixed- 
force-level-breakpoint  battle  with  square-law  attrition. 

The  following  assumptions  yield  the  stochastic  square-law  attrition 
process. 

Al.)  The  attrition  process  depends  only  on  the  current  system 
state  and  time  and  not  on  the  past  history  (this  assumption  is 
usually  referred  to  as  markovian  property). 

«•>  Theprob.b.,.tyP(;“x1«^fV!ioH;Vhe)  ’ 

A3.,  The  Prcbabimy  pfe )  *  » 

A4. )  The  probability  of  more  than  one  casualty  occurring  in 
the  time  interval  t  to  t+h  is  of  the  order  of  magnitude  o(h), 
where  lim^o(h)/h  =  0. 

A5.)  No  more  casualties  can  occur  once  m  =  m^  or  n  =  nbp. 

Making  the  time  interval  h  infinitesimally  small,  the  following 
set  of  forward  CHAPMAN- KOLMOGOROV  equations  can  be  developed.  Let 
P(t,m,n)  be  the  probability  that  the  system  is  in  state  (m,n)  at  a 
time  t.  For  convenience  each  equation  is  related  to  a  region  in  the 
state  space  shown  in  Fig.  2. 
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For  m=mg  and  n=ng,  Region  I 
dP 

^ ^ ( t , nig > n Q )  -  - ( ang+bmg )  P(tjirig,ng)  (2.12) 

for  m^p<m<mQ  and  n  =  ng ,  Region  II 
d  P 

^■(t,m,n0)  =  anQ  P ( t ,m+l ,nQ ) - ( anQ+bm)  P(t,m,n0)  (2.13) 

for  nbP<n<ng  and  fn=mQ,  Region  III 

rl  P 

^(t,m0,n)  =  bmQ  P ( t ,mQ  ,n+l ) - ( an+bmQ )  P(t,mQn)  (2.14) 
for  m.<m<mQ  and  nbp<n<no»  Region  IV 

d  P 

jjr(t,m,n)  =  anP(t,m+l ,n)+bmP(t,m,n+l )-(an+bm)P(t,m,n) 

(2.15) 

for  rc=mbp  and  nbp<n<nQ»  Region  V 


A  p 

dT< 1 ’mbp ,n )  =  anP(t,mbp+l  ,n)  (2.16) 

for  n=n.  and  rnbp<m<n1o »  Region  VI 

^■(t,m,nbp)  =  bmP(t,m,nbp+l )  (2.17) 

for  m=rnbp  and  n=nbp>  Region  VII 

P(t,mbp,nbp)  *  0  for  all  t  (2.18) 

because  of  the  definition  of  a  breakpoint  force  level.  The 
initial  condition  is 

P(0,m0,n0)  =  1.0  (2.19) 
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As  P(t,m,n)  is  a  joint  probability  distribution,  the  follow 
ing  must  also  be  true 

0*  P(t,m,n)-1 .0  (2.20) 

and 


P(t  ,m,n)  =  1 .0  (2.21 ) 

m3mbp  n=nbp 
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III.  SOLVING  THE  DIFFERENTIAL  EQUATIONS 

A.  SOLUTIONS  FOR  THE  DETERMINISTIC  MODEL 

Force  levels  as  a  function  of  time  were  already  given  in  equations 
(2.5)  and  (2.6).  It  is  relatively  easy  to  obtain  analytical  solutions 
for  such  simple  deterministic  models  as  described  before.  On  the  other 
hand,  it  is  necessary  to  point  out  that  for  models  with  any  degree  of 
operational  realism,  analytical  methods  for  solving  the  differential 
equations  are  usually  not  available.  Therefore,  some  numerical  method 
with  a  digital  computer  is  usually  used.  TAYLOR  [10]  has  summarized  in 
his  Appendix  C  the  most  widely  used  numerical  methods,  a  discussion  of 
which  seems  unnecessary  at  this  point. 

B.  SOLUTIONS  TO  THE  STOCHASTIC  MODEL 

Even  for  this  relatively  simple  stochastic  model  with  fixed  force 
level  breakpoints,  which  are  usually  nonzero,  a  complete  set  of  general 
analytical  solutions  for  the  differential  equations  (2.12)  through  (2.19) 
has  not  been  found.  On  the  other  side,  given  some  minor  restrictions 
like  a  fight  to  the  finish  or  equal  attrition  rate  coefficients,  solu¬ 
tions,  or  at  least  solution  methods  have  been  proposed  which  will  be 
briefly  discussed  in  the  next  section. 

First,  the  method  for  getting  the  state  probabilities  used  here  will 
be  described.  Numerical  solutions  were  obtained  using  the  fourth-order 
RUNGE-KUTTA  method,  which  is  probably  one  of  the  best  known  finite  dif¬ 
ference  approximations  to  ordinary  differential  equations  (next  to  the 
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EULER-CAUCHY-method) .  To  increase  the  accuracy  of  the  overall  solutions 
analytical  results  for  certain  regions  of  the  state  space  were  substituted. 
These  analytical  solutions  will  be  stated  now.  For  region  I,  i.e.  no 
casualties  on  either  side,  the  solution  to  (2.12)  is 

P(t,m0,n0)  =  exp  -(an0+bmQ)t  (3.1) 

which  can  be  derived  by  the  standard  method  used  for  this  kind  of  dif¬ 
ferential  equation.  For  the  boundary  cases,  region  II  and  III,  i.e. 
one  of  the  two  sides  has  not  yet  had  a  casualty,  TAYLOR  [10]  has  given 
the  analytical  expression  as 
for  mbp<m<m0  and  n=nQ,  Region  II 

P(t,m,nQ)  =  j,  janQ/b(ebt-l  )| J  exp(- (bm0+anQ)  t)  ,  (3.2) 

where  J  =  irig-m 

for  nbp<n<nQ  and  ,  Region  III 

P(t,mQ,n)  =  |bmQ/a(eat-l ) |K  exp ( - ( bmQ+an0 )  t)  ,  (3.3) 

where  K  =  nQ-n. 

These  two  equations  (3.2)  and  (3.3)  were  obtained  by 
recursively  solving  equations  (2.12),  (2.13)  and  (2.14)  "from 
the  top  down." 

C.  OTHER  ANALYTICAL  SOLUTIONS 

The  solutions  or  solution  methods  for  getting  the  state 
probabilities  will  only  be  stated  for  the  square-law  attri¬ 
tion  process.  Only  two  were  used  for  the  numerical  work  for 
this  thesis. 
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Apparently,  one  of  the  "oldest"  analytic  solutions  was 
given  by  BROWN  [3]  in  1955  for  the  general  stochastic 
LANCHESTER-type  combat  model  with  time  independent  attrition 
rates.  His  approach  and  solution  will  be  briefly  outlined 
for  square-law  attrition.  Consider  a  path  from  state 
(  mQ  ,  n0  )  to  some  state  (m,n).  This  path  can  be  described 
as  a  sequence  of  J  =  (  mg  -m)  zeros  and  K  =  (  ng  -n)  ones, 
where  a  zero  corresponds  to  an  X  casualty  and  a  one  to  a  Y 
casualty.  Using  the  binary  representation  of  a  positive 
integer,  one  can  define  to  each  realization  of  a  battle  path 
an  integer  k  given  by 


k  =  dk  ,  1 dk  ,2 ‘ ’ ’ dk  ,  J+K 

where  dk  r  *  1  if  the  rth  casualty  along  a  battle  path 
corresponding  to  k  is  a  Y  casualty  and  dk  r  *  0  otherwise. 
Also  let  Ij  K  be  the  set  of  all  positive  integers  whose 
binary  representation  contains  exactly  K  ones  and  J  zeros. 
Then 


"k,r  ‘  Vr* 


t*ki 

j=l  K’J 


k,r 


r 


Then  BROWN  [3]  has  shown  that 
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i  _  2  J  +  K-l 

n)  =  T7  21  _/  IT-  Kk, 

keIJ,K  ""  r=0 


exp(-iut)-exp(-tl (mk  >r ,nk  >r) ) 
l{mk,r’nk,rMu 


du  ,  (3.4) 


where  i  =  if- 7  , 

l(m,n)  «  an+bm, 

K  .  9k ,r+1 _ 

k,r  ~  1  -  i  u  /  1(m|c>r»nk  >r) 

and 

gk,r+l  =  dk,r+lank,r+(1'dk,r+l )bmk,r  * 

There  was  no  indication  that  nonzero  breakpoints  were 
excluded.  A  discussion  of  this  solution  follows  in  the  next 
section . 

About  14  years  later,  in  1969,  CLARK  [4]  proposed  another 
approach  which  TAYLOR  [10]  called  a  "hybrid  analyti cal -numer¬ 
ical"  method.  The  restriction  is  that  the  breakpoints  have 
to  be  zero,  i.e.  it  is  a  fight  to  the  finish.  Although 
proposed  for  a  general  time  independent  attrition  function, 
this  approach  will  be  outlined  for  "square-law"  attrition. 
Then  according  to  CLARK  [4]  the  state  probabilities  are  given 
by 
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for  0<mSm0  and  Qcning 


m  n  n 


P(t,m,n)  =  ‘exP  (-(ak+bj )tj  ,  (3.5) 

j=m  k=n 


for  0<m<m»  and 


n  =  0 


mo  no 


P(t,m,0)  =  cj)’°  +  Z.  21  C;*°  •exp(-( ak+bj )t)  ,  (3.6) 

u,u  j=m  k  =  l  J,k  v  ' 


for  m=0  and  0<n<n, 


m0  nQ 


P( t ,0 ,n )  =  C°>2  +  Z  Z.  •  exp  (- ( ak+b  j  )  t)  ,(3.7) 


and  at  last  there  is  to  remember  that  P(t,0,0)  =  0  for  all 
times  t.  The  constants  are  determined  by  a  system  of 

J  j  K 

partial  difference  equations. 


For  0<m<j<m0  and  0<n<k<nQ 


rm,n 
j » k 


rm+l ,n  ,  rm ,n+l 
an  ,k  bm  ,k 

a(n-k)+b(m- j) 


9 


(3.8) 


I 


Also  for  0<mSjsmQ  and  0=n<k£ng 


,m  ,0 


ak+b  j 


(3.12) 


similarly  for  0=m<j<mg  and  0<n<k<ng 


r0,n 
Lj,k  = 


an  J,k 
ak+bj 


(3.13) 


Then  for  l<m<mr 


mQ  n. 


j-m  k=l 


(3.14) 


and  finally  for  iSnSng 


v  u 

,n  _  ^T“  V"  c0,n 

,0  Z-  Z_  Uj,k 

j=l  k=n 


(3.15) 


Though  having  the  publishing  date  of  1979,  the  next 
approach  was  published  in  June  1980  by  8ILLARD  [1J.  She 
considered  the  LANCHESTER-type  square-law  attrition  combat 
model  as  a  pure  death-process  and  applied  SEVERO's  [8] 
recursive  theorem  for  solving  differential  equations.  As 
before,  only  a  fight  to  the  finish  has  been  considered. 
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The  first  step  is  to  identify  each  point  (m,n)  in  the 
state  space  by  a  counting  coordinate  k,  where 
k  =  (  mQ+l ) (  nQ+l ) 

-m(  nQ+l)-n  .  (3.16) 

Then 

P [t,m,n] =  P[t,k]  (3.17) 

and  the  differential  equations  (2.12)  through  (2.18)  take  on 
a  slightly  different  form.  As  an  example,  (2.15)  will  be 
given  by 

^r(t,k)  =  anP [t ,k-nQ-l ] +bmP  [t  ,k-l ] 

-(an+bm)P  [t  ,k]  .  (3.18) 

The  whole  set  of  differential  equations  was  then 
expressed  in  matrix  terms  as 

afp(t)  =  BP(t)  (3.19) 

with  a  solution  given  as 

P(t)  =  Ce  (t)  (3.20) 

where  e^(t)  is  the  (mQ+l  )  ( nQ+l ) xl  -vector  with  elements 
exp(bkt)  with  bk  being  the  kth  diagonal  element  of  the 
matrix  B.  The  matrix  can  be  partitioned  into  submatrices, 
whose  m-coordinate  is  common,  due  to  the  ordering  defined  by 
the  counting  coordinate  k  (Equation  3.16). 

Then 

B.  =  ( b^  y)’  u  ,v  =  1,2..  .mg+1 
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where  the  submatrices  buv  have  the  elements 
— 'U  V  ""  P  ,  q  ~  1,2...  Oq+1  . 

So  the  matrix  B_  has  the  elements 

buu(p,p)  =  ~  a ( n0" P+1 )-b(m0-u+l ) 

for  u  =  1,2. ..mg  and  p  =  l,2...ng 

buu(p,P-D  =  b(mQ-u+l)  for  u  =  l,2...m0+l 

p  =  2 , 3 . . . nn+l 

and  u 

bu  ,u-l (p,p)  =  a(nQ-p+l)  for  u  =  2,3...m0  +  l 

P  =  1,2...  r>g  +  l  . 

All  other  elements  are  zero. 

Thus,  the  matrix  B_  has  at  most  3  nonzero  entries  per  row 
or  column.  The  matrix  £  can  be  partitioned  in  the  same  way. 
Then  using  SEVERO’s  [8]  theorem  and  the  special  form  of  the 
matrix  B^,  only  a  part  of  the  C_-matrix  needs  to  be  determined. 
This  part  will  be  omitted  here,  but  the  final  result  will 
be  given  by 

k 

P(t,k)  =  H  c(k,j) .exp(bjt)  (3.21) 

j  =  l 

X.  L. 

where  c(k,j)  is  the  (k,j)  element  of  the  solution  matrix  C. 

The  previous  two  approaches  have  required  that  the  force 
level  breakpoints  be  zero.  Now,  a  result  will  be  given  whose 
restriction  is  that  the  attrition  rate  coefficients  be  equal, 
but  nonzero  breakpoints  are  allowed.  For  further  reference 
it  will  be  called  the  Equal -Attri ti on-Rate-Coef f i ci ent- 
Solution  (EARCS). 
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Let  a  =  b  =  f. 

For  mbp<m<mQ  and  nbp<n<nQ 


PCt’m-n)  =  Tvi^-hTl 


( 1-e  )  0 


ft^n+n0-m-n 


exp(-f (m+n)t) 


(3.22) 


for  m  =  mbp  and  nbp<n<nQ 


P ( t  *  mb  p  * n )  = 


k=0 


fn  C(mbp-H,n) 


(!)  P 


J! 

-exp(-ft(mbp 
f (mbp+l +n+k ) 


+1+n+k 


)i 


}■ 


(2.23) 


where  J  =  mo+no"mbp“^ “n 
f°r  "  =  nbp  and  mbp<m-mQ 


P(tfm,nb  )  = 


fm  C(m,n  bp+1> 
K! 


.  /  k  \  J"1-exp(-ft(m+Unbp+j))’t 

Z-(-1)J  \j  )  l  f  (m+1  +nbp+j  )  ^  J’ 


(3.24) 


where  K  =  mQ+nQ-m-nbp-l  . 

The  coefficients  C(m,n)  satisfy  for  ^bp<m<t”g  and  nbp<n<n0 
the  following  partial  difference  equation 

C(m,n)  -  nC(m+l,n)  +  mC(m,n+l)  (3.25) 
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with  the  boundary  conditions 

mn-m 

C(m,ng)  =  (ng)  for  < m<m 


and 


nQ-n 

C(m0,n)  =  (mg)  for  nbp<n<nQ 

This  result  has  been  developed  using  a  method  verbally  pro¬ 
posed  by  TAYLOR.  The  method  will  now  be  outlined.  ISBELL 
and  MARLOW  [6]  described  a  stochastic  LANCHESTER-type  attri¬ 
tion  process  with  a  different  attrition  function.  Instead  of 
the  attrition  of  one  force  being  proportional  to  the  number  of 
of  enemies  of  the  opposing  force  as  in  the  square-law  attri¬ 
tion  (e.g.  for  the  attrition  of  the  M-force  let 
A(m,n)  =  an 

be  the  attrition  rate  and 
B(m,n)  =  bm 

the  attrition  rate  for  the  N-force  with  square-law  attrition), 
their  attrition  rates  looked  like 
A(m,n)  =  an+cm 


and 

B(m,n)  -  bm+dn 
with  the  restriction  that 
a+c  3  b+d. 

But  with  c=d=0  and  a=b=f  we  are  back  to  square-law  attrition. 
This  leads  to  equations  (3.22)  and  (3.25).  Equations  (3.23) 
and  (3.24)  were  derived  in  the  following  way  (e.g.  for  3.23). 
Solving  equation  (3.22)  for  m  =  mbp+l  yields 
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1 


v 


P(t,mbp,n)  = 


C (mbp+1 *n> 

tmO+nO-rabp’1‘n,! 


exp 


|-ft(m 


bp+1+n>) 


.  (i.e-ft)mO+nO‘tT'bp'1‘n 


(3.26) 


substituting  for  the  second  factor  its  BINOMIAL  expansion 

0  / 1  \ 

,k 


e-ftk(1)J-k 


k=0 


with  J  *  n,o+nO”mbp"^  anc*  mu^tiPlyin9  through  by  the  third 
factor.  Then  using  equation  (2.16),  the  differential  equation 
for  m  =  mbp  and  nbp<n-no  »  and  substituting  equation  (3.26) 
into  the  extended  form,  it  can  now  easily  be  integrated  to 
yield  equation  (3.23). 


D.  DISCUSSION 

In  the  discussion  of  the  analytical  solutions  outlined  in 
the  last  section,  there  is  one  important  point.  BROWN  [3] 
himself  points  out  that  unless  m  is  close  to  mQ  and  n  is 
close  to  nQ  his  result  (equation  (3.4))  is  of  "little  practi¬ 
cal  interest."  Most  of  the  analytical  solutions ,  especial ly  for 
more  general  LANCHESTER-type  models,  have  little  more  than 
"symbolic"  character.  BROWN'S  solution  is  a  good  example 
of  that. 

In  comparing  the  solutions  given  by  BILLARD  [2]  and 
CLARK  [4],  this  author  has  the  feeling  that  both  solutions  are 
equivalent  and  only  the  representation  is  different.  This 
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intuitive  guess  needs  verification.  It  may  be  concidence 
that  CLARK  [4]  and  SEVERO  [8]  published  their  work  in  the 
same  year. 

The  last  presented  solution  (equation  (3.22)  through 
(3.24))  seems  to  be  relatively  handy  for  use  on  a  computer. 

It  has  a  big  advantage  over  numerical  solution  methods  other 
than  its  accuracy,  because  it  is  an  exact  result.  Like 
CLARK's  method,  the  coefficients  have  to  be  calculated  only 
once  for  a  given  set  of  input  data.  Then,  to  get  the  state 
probabilities  for  a  certain  point  in  time  you  have  to  make 
only  one  set  of  calculations,  as  opposed  to  the  numerical 
methods  where  one  has  to  go  from  time  t  =  0  to  time  t  =  t 
in  small  time  steps  and  then  have  only  an  approximate  result. 

E.  IMPORTANCE  OF  THE  STATE  PROBABILITIES 

The  state  probabilities  as  a  function  of  time  are  the 
key  to  calculating  several  quantities  of  interest.  These 
are  expected  force  levels  as  a  function  of  time,  variances 
and  standard  deviations  in  the  force  levels  and  also  the 
probability  of  winning.  These  quantities  are  necessary  to 
legitimately  compare  the  stochastic  with  the  deterministic 
results . 

To  get  at  least  a  feeling  of  how  the  state  probabilities 
evolve  over  time,  the  joint  probability  distribution  will  be 
presented  in  a  3-D-picture.  It  is  indeed  surprising  that 
more  use  has  not  been  made  of  computer  graphics  to 


investigate  the  dynamics  of  a  stochastic  LANCHESTER-type 
combat  model.  Table  1  gives  the  data  used  for  the  next  five 
figures.  These  figures  may  be  thought  of  as  "snapshots"  of 
the  joint  probability  for  the  survivors  in  this  battle  taken 
at  a  sequence  of  increasing  times. 
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TABLE  1 


Data  for  the  Numerical  Example  1 


Force  Levels 

3 

O 

II 

4* 

O 

o 

II 

Q. 

-Q 

E 

o 

*** 

II 

o 

c 

"bp  ■  ° 

Attrition  Rates 

a  =  0.008 

M  casualties 
minute  and  N 

per 
f  i  rer 

b  =  0.004 

N  casualties 
minute  and  M 

per 
f  i  rer 

At  Times 

t1  =  0.025 

tf 

t2  =  0.250 

tf 

t3  =  0.500 

tf 

t4  =  0.750 

tf 

t5  =  1.000 

tf 

where  =  155.81  minutes  is  the  time  a 
deterministic  battle  with  the  same  force 
levels  and  unit  breakpoints  ends,  i.e. 
x(tf)  =  0.0  and  y(t^)  =  28.28 
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Figura  3 


-  FLOT  OF  JOINT  PROBABILITIES  P(t,a,n; 


with  data  according  to  Tabla  1 


rigur*  4 


-  rLOT  OF  JOINT  PROBABILITIES  P(t,a 


with  data  according  to  xabla  1 


Figure  5 


-  .PLOT  OF  JOINT  PROBABILITIES  P(t 


with  data  according  to  Tabla  1 


rigurw  6  -  PLOT  OF  JOINT  PROBABILITIES  P(t,m,n) 
with  data  according  to  Tablw  1 
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F(t,m,n) 


rigura  7  -  PLOT  OF  JOINT  PROBABILITIES  P(t,a,n) 
with  data  according  to  Tabla  1 
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At  the  time  t  =  0  all  probability  is  located  at 
(m0,n0)  in  the  state  space  (Region  I).  As  time  passes,  the 
probability  mass  is  distributed  over  more  states,  with  the 
mode  moving  away  from  the  diagonal  towards  the  winning  side. 
All  points  in  the  state  space  with  breakpoints,  i.e.  ( p , n ) 
and  ("1»n^p)  for  all  m  and  n,  are  absorbing  states,  probability 
mass  is  absorbed  in  that  states.  The  sum  of  probability  mass 
in  Region  V  (see  Fig.  2)  represents  the  probability,  that  the 
N-force  wins  at  that  given  time,  in  Region  VI  that  the  M-force 
wi  ns . 

The  next  five  figures  show  a  similar  sequence  of  plo„s  for 
the  joint  probability  P(t,m,n)  for  the  force  levels  M(t)  and 
N(t).  The  data  is  explained  in  Table  2.  Note  the  small  dif¬ 
ferences  because  of  the  nonzero  force  level  breakpoints. 
Probability  mass  having  reached  the  breakpoint  "piles"  up 
there.  The  state  space  is  reduced  by  the  fixed  force  level 
breakpoints,  but  the  probability  distribution  evolves  basical¬ 
ly  in  the  same  qualitative  manner  as  in  the  previous  example. 
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1 

1 

TABLE  2 


Data  for 

the 

Numerical  Example  2 

Force  Levels 

3 

O 

II 

-F* 

O 

m.  =  8 
bp 

o 

II 

o 

c 

nKr.  =  8 
bp 

Attrition  Rates 

a  =  0.008 

M  casualties 
minute  and  N 

per 
f  i  rer 

b  =  0.004 

N  casualties 
minute  and  M 

per 
f  i  rer 

At  Times 

t1  =  0.025 

tf 

t2  =  0.250 

tf 

t3  =  0.500 

tf 

t4  =  0.750 

tf 

t5  »  1.000 

lf 

where  t^  = 

120 

.68  minutes  is 

the  time  a 

deterministic  battle  with  the  same  force 
levels  and  unit  breakpoints  ends,  i.e. 
x(tf)  =  8.0  and  y(t^)  =  28.83 
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P(t,B,n) 


Figure  8  -  t'LOT  OF  JOINT  PROBABILITIES  P(t,a,n) 

with  data  according  to  Table  2 
19 


EM 


figura  1u  -  PLOT  OF  JOINT  PROBABILITIES  P(t,a,n; 
with  data  according  to  labia  2 


f'igura  11  -  PLOT  OF  JOINT  PROBABILITIES  P(t,a, 
with  data  according  to  Tabla  2 
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rigur*  12  -  PLOT  OF  JOIKT  PROBABILITIES  P(t,a,n) 
vith  data  according  to  Tabla  2 


IV.  EVOLUTION  OF  THE  FORCE  LEVELS 


A.  THE  DETERMINISTIC  MODEL 

The  evolution  of  force  levels  as  a  function  of  time,  X(t)  and  Y(t), 
has  already  been  stated  in  equations  (2.5)  and  (2.6).  The  next  two 
figures  show  the  force  levels  with  different  breakpoints.  It  is  easy  to 
realize  that  introducing  a  nonzero  breakpoint  does  not  change  the  under¬ 
lying  function.  In  other  words,  using  X(t)  as  an  example,  for  both 
figures  the  same  curve  was  used  but  at  the  point  where  the  X-force  reaches 
its  breakpoint,  the  curve  is  "cut.”  From  that  point  in  time,  there  are 
no  more  changes  in  the  force  levels.  So  introducing  a  nonzero  break¬ 
point  only  shifts  the  discontinuity  (marked  by  DX1  in  Fig.  13)  up  along 
the  curve  to  DX2. 

The  probability  for  one  side  to  win  is  either  one  or  zero,  because 
it  is  a  deterministic  model.  To  easily  determine  which  side  is  going  to 
win,  a  victory  prediction  condition  can  be  obtained  by  solving  each  force 
level  equation  (2.5)  and  (2.6)  for  the  time  to  reach  its  breakpoint  txbp 
by  substituting  xU)=xbp  and  tybp  by  substituting  Y(t)=ybp.  Then  X  will 
win  if  tybp  <  txbp,  which  leads  to  the  prediction  condition.  X  will  win 
a  fixed  force  breakpoint  battle  if  and  only  if 


This  shows  that  given  the  initial  data  one  can  predict  the  outcome 
of  the  battle  in  terms  of  force  levels  and  time  until  the  battle 
finishes. 
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tirne  t 


Figure  14  -  rORCE  LEVELS  OVER  TIME 


Deterainistic  Model  with  Breakpoint  x.  *  0 


B.  THE  STOCHASTIC  MODEL 


Every  possible  state  (m,n)  in  the  state  space  has  associated  with  it 
a  certain  probability  between  zero  and  one  which  is  a  function  of  time. 

In  order  to  gain  more  insight  into  the  stochastic  process  not  a  single 
realization  of  a  battle  has  to  be  considered,  but  an  average  battle. 
Therefore  the  expected  value  of  the  force  levels  (i.e.  averages)  as  a 
function  of  time  and  the  variances  in  the  force  levels  were  investigated. 

The  straight  forward  way  to  compute  the  expected  force  levels  involves 
the  knowledge  of  the  probability  distribution  P(t,m,n).  Then 

mQ  n0 

E(M(t)  )=  P(t,m,n)  (4.2) 

mbp  nbp 


n0 

E(N(t)  )=  ^  n 

nbp  mbp 

are  the  expected  values  of  the  force  levels  as  a  function  of  time. 

There  are  some  other  ways  to  calculate  the  expected  force  levels, 
one  of  which  will  be  stated  here.  Recall  the  "hybrid-analytical- 
numerical"  method  to  get  the  state  probabilities  (equations  3.5  through 
3.15).  CLARK  [4]  has  also  shown  that  the  itfl  moment  of,  for  example, 
the  M-force  level  may  be  computed  as 


t 


P(t,m,n) 


(4.3) 


E(Mi(t) 


)  • 


exp(-(ak+bj)t  ) 


(4.4) 


with  for  1-j-mg  and  ISk^ng 
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(4.5) 


3(i)  = 

Jj,k 


m 


m=l 


m,n 

j,k 


r0,n 

^0,0 


(4.6) 


This  again  emphasizes  the  strong  point  of  CLARK's  solution.  For  a 
given  set  of  battle  parameters  the  coefficients  cl?’!)  have  to  be  computed 
only  once.  Then  with  this  information  and  relatively  small  computa¬ 
tional  effort  not  only  the  state  probabilities  but  also  the  first  and 
second  moment  of  the  force  levels  can  be  computed.  This  determines  the 
variance  in  the  force  levels  at  the  same  time,  e.g. 

Var(M(t)  )  =  E(M2  (t)  )  -  E(M(t))  •  E(M(t))  .  (4.7) 

On  the  other  side,  the  weak  point  is  that  CLARK  [4]  considered  only 
breakpoints  equal  zero. 

Many  authors  have  discussed  one  side's  probability  of  winning  or 
probability  of  winning  conditioned  on  a  certain  number  of  survivors, 
which  always  eliminated  the  parameter  time  by  integrating  from  time  t=0 
to  infinity.  This  may  be  legitimate  to  answer  absolute  (meaning  time 
independent)  questions  about  who  will  win,  but  for  direct  comparisons 
with  the  deterministic  model,  this  author  has  the  feeling  that  the  best 
question  to  ask  regarding  a  winner  is: 

What  is  the  probability  of  one  side  winning  given  the  stochastic  battle 
lasted  as  long  as  the  deterministic  one? 

The  calculation  of  these  probabilities  gives  another  interesting 
probability,  because  given  the  time  t  =  t^  (final  time  of  the  determin¬ 
istic  battle) 
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P(battle  has  not  yet  finished)  = 

1-P(M  winsj t=tf)-P(N  wins j t=t^)  .  (4.8) 

There  has  also  been  work  done  regarding  the  distribution  of  the 
time  to  finish  a  battle.  But  this  is  beyond  the  scope  of  this  thesis 
(SPRINGALL  [9)). 

C.  DIFFERENCES  IN  THE  FORCE  LEVELS 

The  deterministic  model,  especially  the  force  level  equations  (2.5) 

and  (2.6),  describe  a  process  with  continuous  state  parameters  where,  in 

reality,  the  possible  states  are  integers.  Quoting  from  LANCHESTER  [7]: 

Since  th  •  forces  actually  consist  of  a  finite  number  of  finite  units 
(instead  of  an  infinite  number  of  infinitesimal  units)  the  end  of  the 
curve  must  show  discontinuity,  and  break  off  abruptly  when  the  last 
man  is  reached;  the  law  based  on  averages  evidently  does  not  hold 
rigidly  when  the  numbers  become  small. 

LANCHESTER  suggested  that  his  differential  equations  may  be  good 
approximations  only  as  long  as  the  force  sizes  are  large.  He  also  stated 
that  the  equations  are  based  on  averages,  implying  an  underlying  stochas¬ 
tic  process. 

This  shows  that  there  must  be  a  difference  in  the  force  levels  which 
should  become  significant  when  the  number  of  combatants  is  small.  This 
difference  was  called  bias  by  CLARK  [4]  and  TAYLOR  [10].  It  can  be  shown 
that 

-aE(N(t))+aBn(t)  (4.9) 

and 

d|CaftI}-  ■  -bE(M(t))+HBra(t) 


where 


(4.10) 


t»mbp,n) 


(4.11) 


i 


Bn(t)  s  "bp  2_  P(t-m-nbp)+  2_  nP( 
"bp*’  "bp*' 


and 


Bm(t>  =  mbp2_  P(t,mbp,n)+  2_  mP(t,m,nbp)  .  (4.12) 


V1 


V1 


The  bias  terms  Bm(t)  and  Bn(t)  can  be  interpreted  as  the  expected 
values  of  M(t)  or  N(t)  conditioned  on  the  fact  that  the  battle  has 
already  ended  at  time  t,  for  example, 

Bm(t)=E  [M(t))|N(t)=nbp  or  M(t)«n^p]  (4.13) 

In  other  words,  equation  (4.9)  says  the  expected  casualty  rate  of  the 
M- force  is  proportional  to  the  expected  number  of  survivors  of  the 
N-force  given  neither  of  the  two  forces  has  reached  its  breakpoint. 

Define  the  bias  of  the  X-force  as  Ax(t)  =  E(M(t))-X(t)  and  the  bias 
of  the  Y-force  as  Ay(t)  =  E(N(t))-Y(t) .  Then  using  equations  (2.8)  and 
(2.9)  together  with  (4.9)  and  (4.10)  it  follows  that 

Ax  =  -aAy+aBn(t)  (4.14) 

^  Ay  =  -bAx+bBm(t)  (4.15) 

with  the  initial  conditions  x(0)=)  and  Y(0)=).  This  has  the  solution 

t 

Ax(t)=  VS  /  IBn(s) */a/bcosh (t-s) 

-Bm(s)sinh( ~\fab  ( t-s ) ) >ds ,  (4.16) 
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and 


t 

Ay(t)=  Vab  J  (Bm(s)  Vb/acosh(  Vab(t-s) 

-Bn(s)sinh( Vab(t-s))}ds  .  (4.17) 

Since  for  a  fixed,  nonegative  argument  z,  the  cosh(z)  is  always 
greater  than  the  sinh(z),  it  is  easy  to  visualize  that  in  most  of  the 
cases  both  biases  are  positive,  meaning  the  expected  force  levels  of  the 
stochastic  model  are  higher  than  the  deterministic  force  levels.  This 
has  been  shown  by  CLARK  [4]  and  CRAIG  [5)  and  confirmed  by  this  author. 
In  the  rest  of  the  cases  the  winner's  bias  is  negative  or  close  to  zero 
and  the  loser's  bias  is  positive.  Two  examples  are  given  in  Table  3. 
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TABLE  3 


Examples  for  cases  where  one  bias  is  positive  and  the  other 
bias  is  negative. 


A.  Y-force  wins  in  a  fight  to  the  finish 


xo  =  y0  =  40 

xbP  *  0 

^bp  '  ^ 
tf  =  15.581 


a  =  0.08 
b  =  0.04 


Ax(tf)  =  3.22 
Ay(tf)  =  -0.23 


B.  Y-force  wins  in  a  fight  with 
but  different  breakpoints. 

xQ  =  y0  =  15 


x 


bp 


12 


y 


bp 


6 


tf  =  2.85 


equal  initial  force  levels, 

a  =  0.08 
b  =  0.08 

Ax ( tf)  =  0.15 
Ay(t^)  =  -0.10 
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An  interesting  point  has  to  be  mentioned  regarding  case  B  of 
Table  3.  Here  the  expected  force  level  of  the  winner  is 
smaller  than  the  expected  force  level  of  the  loser  at  the 
time  a  deterministic  battle  ends. 

To  get  a  better  feeling  for  the  differences  in  the  force 
levels,  the  next  two  figures  show  as  an  example  a  large  spec¬ 
trum  of  force  level  behavior.  The  data  and  notation  is 
described  in  Table  4.  There  are  four  battles  with  four  dif¬ 
ferent  breakpoints  drawn  as  they  evolve  over  time  from  zero 
to  the  time  an  equivalent  deterministic  battle  ends. 
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TABLE  4 


xo 

=  40 

V> 

=  40-0. 2i 

Yo 

=  40 

V'> 

II 

-P* 

O 

1 

o 

ro 

a 

=  0.09 

for  i 

=  1,2,3, 4 

b 

=  0.07 

Mi 

=  E  M(t) 

for  i 

=  1  ,2,3,4 

Ni 

=  E  N  ( t ) 

Where  the  index  i  corresponds  to  the  battle  with 

X.  L. 

the  il  breakpoint  force  level: 

X-force  is  always  loser 
V-force  is  always  winner 
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CRAIG  [5]  has  formulated  hypotheses  concerning  the  biases  in  the  average 
force  levels  based  on  his  work.  His  hypotheses  were  partially  confirmed, 
but  in  some  cases  they  have  to  be  modified.  Therefore,  another  similar 
set  of  hypotheses  will  be  given  and  supported  by  Fig.  15  and  Fig.  16,  as 
well  as  some  of  the  later  figures. 

H  1)  Given  fixed  initial  force  levels  and  attrition  rate  coefficients, 
as  the  breakpoint  force  levels  increase,  the  numerical  bias  for  the 
loser  decreases.  The  biases  for  the  winner  do  not  show  this  monotone 
trend  except  for  the  case  of  symmetric  parity. 

H  2)  Everything  else  constant,  the  bias  of  the  loser  increases  with 
increasing  initial  force  levels;  this  is  also  true  in  the  symmetric 
parity  case  for  both  forces. 

H  3)  Given  the  initial  force  level  ratio  is  close  to  one  at  the  time 
corresponding  to  the  end  of  the  deterministic  battle,  the  bias  of  the 
loser  is  always  larger  than  the  bias  of  the  winner. 

H  4)  At  the  time  corresponding  to  the  end  of  a  deterministic  battle, 
the  biases  become  larger  as  the  forces  come  closer  to  parity. 

H5)  The  biases  at  times  corresponding  to  less  than  one  half  the 
duration  of  the  deterministic  battle  are  negligible. 

The  case  of  symmetric  parity,  i.e.  equal  initial  force  levels, 
breakpoints  and  equal  attrition  rate  coefficients  seems  to  be  kind  of  a 
"limiting"  case.  For  example,  at  parity  the  biases  of  both  forces  behave 
in  the  same  manner  and  are  equal.  In  Fig.  17  the  biases  at  the  time  a 
deterministic  battle  ends  as  a  function  of  the  initial  force  levels  and 
as  a  function  of  the  breakpoints  are  presented.  It  is  also  another 
verification  for  the  hypotheses  H  1  and  H  2. 
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Fi*ur*  18  -  BIASES  WITH  DIFFERENT  INITIAL  FORCE  LEVELS 


Figur#  19  -  BIASES  WITH  DIFFERENT  ATTRITION  RATE 


COEFFICIENTS 


Considering  the  changes  in  magnitude  of  the  biases  for  battles  like 
in  Fig.  15,  16  and  17,  CRAIG  [5]  came  to  the  conclusion  that  when  the 
forces  are  closer  to  parity,  the  biases  at  the  deterministic  battle's 
end  increase.  Several  battles  were  fought  starting  with  symmetric  parity 
and  then  varying  the  force  levels  and  the  attrition  rate  coefficients. 
Sample  results  are  shown  in  Fig.  18  and  Fig.  19.  In  Fig.  18,  the  initial 
force  levels  were  changed  giving  the  Y-force  ten  and  twenty  percent 
higher  initial  force  levels.  The  biases  at  the  end  of  the  deterministic 
battle  are  plotted  as  a  function  of  the  breakpoint  force  level  ratio 

r  =  xbp/x0* 

Fig.  18  supports  CRAIG's  hypothesis,  as  does  Fig.  19.  Here  not  the 
initial  force  levels  but  the  attrition  rate  coefficients  were  changed  in 
order  to  deviate  from  symmetric  parity.  The  last  way  to  deviate  from 
symmetric  parity  is  a  case  where  the  deterministic  model  gives  equal 
answers  for  different  battles.  The  data  and  the  results  are  shown  in 
Table  5. 
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TABLE  5 


Equal  initial  force  level  battle  with  non  equal 
XQ  =  15  Y0  =  1 5  a  *  b  =  f  = 


X, 

bp 

Ybp 

Ax(tf 

12 

12 

1.08 

12 

9 

0.71 

12 

6 

0.15 

12 

3 

0.68 

9 

9 

1.62 

9 

6 

1.14 

9 

3 

1  .06 

6 

6 

2.15 

6 

3 

1.67 

3 

3 

2.84 

breakpoints 

0.08 

Ay(tf ) 


1  .08 
0.63 
-0.10 
0.58 
1  .62 
0.93 
0.75 
2.15 
1.22 
2.84 
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The  X-force  level  at  time  t^  is  the  same  as  the  Y-force  level.  Thus, 
the  expected  force  level  for  the  winner  (i.e.  the  Y-force)  is  smaller  than 
the  expected  force  level  of  the  loser  (the  X-force),  because  Ax(tf)  is 
always  larger  than  Ay(tf).  This  is  easy  to  see  when  one  remembers  the 
way  the  expected  force  level  is  computed  (equation  (4.3)).  Since  the 
Y-force  has  a  lower  breakpoint,  there  are  states  (m,n)  possible  where 
nbp<n<mbp’  Apparently  these  states  have  a  nonzero  probability  associated 
with  them,  which  decreases  the  expected  force  level  below  the  expected 
value  for  the  X-force.  This  might  be  a  starting  point  for  further 
studies. 

D.  VARIABILITY  IN  THE  FORCE  LEVELS 

Naturally  in  the  deterministic  case  there  does  not  exist  any  variabil¬ 
ity  in  the  force  levels.  On  the  other  hand,  for  the  stochastic  model,  the 
variance  in  the  force  levels  as  a  function  of  time  is  a  measure  of  the 
dispersion  of  the  number  of  survivors  about  their  mean  value. 

CLARK  [4]  has  hypothesized  two  different  types  of  behavior  for  the 
variance  in  the  force  levels,  which  are  shown  in  Fig.  20  for  the  data 
presented  in  Table  2.  The  first  type  of  behavior  is  that  of  the  variance 
for  the  N-force,  i.e.  the  variance  increases  monotonely  as  a  function  of 
time  and  is  asymptotic  to  a  limiting  value.  It  was  found  that  this  type 
of  behavior  occurs  when  the  side  is  going  to  win  and  for  the  case  of 
symmetric  parity.  The  second  type  of  behavior  shown  is  the  variance  of 
the  M-force,  Var(M(t)),  as  a  function  of  time.  This  increases  to  a 
maximum  value  then  decreases  asymptotically  to  a  limiting  value.  This 
type  of  behavior  is  associated  with  the  loser  of  the  battle. 
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The  variance  in  the  force  levels  is  a  function  of  the  initial  and 
breakpoint  force  levels,  the  battle  time  and  the  attrition  rate 
coefficients.  Unfortunately  one  does  not  know  what  this  dependence  is. 
Based  on  many  numerical  results  a  set  of  hypotheses  will  be  stated  and 
the  next  figures  will  support  them. 

H  1)  Given  fixed  initial  force  levels  and  attrition  rate  coefficients, 
as  the  breakpoint  force  levels  increase  the  variance  of  the  force  level 
decrease. 

H  2)  Everything  else  constant  the  variance  of  the  loser's  force  level 
increases  with  increasing  initial  force  levels.  This  is  also  true  for 
both  variances  in  the  case  of  symmetric  parity. 

H  3)  Given  the  initial  force  ratio  is  close  to  one  at  the  time 
corresponding  to  the  end  of  the  deterministic  battle  the  variance 
of  the  loser's  force  level  is  always  smaller  than  the  winner's  variance 

H  4)  At  the  time  corresponding  to  the  end  of  the  deterministic  battle 
the  variance  in  the  loser's  force  level  increases  as  the  forces  come 
closer  to  parity.  This  trend  is  not  true  for  the  variance  of  the  win¬ 
ner,  except  for  the  case  of  symmetric  parity. 

It  was  the  intention  of  this  author  that  the  set  of  cases  to  illus¬ 
trate  the  hypotheses  are  the  same  as  in  the  illustrations  (Fig.  17,  18 
and  19)  of  the  hypotheses  about  biases.  So,  Fig.  21  shows  the  force 
level  variances  for  the  different  initial  force  levels  and  different 
breakpoints  for  the  case  of  symmetric  parity.  In  Fig.  22  the  variances 
for  a  battle  with  equal  attrition  rate  coefficients  but  varying  initial 
force  levels  show  that  H  4  is  only  true  for  the  loser.  This  point  is 
emphasized  by  Fig.  23,  where  with  constant  and  equal  initial  force  levels 
the  attrition  rate  coefficients  were  varied. 
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/.0~  to  f=0.0 8 
Case  I  y0=  10 

J£  y0-ii 
III  y0^  1Z 


o.t  aa  Ts  yw>4b 


ERENT  INITIAL 


FORCE  LEVELS 


•  23  -  VARIANCES  WITH  DIFFERENT  ATTRITION  RATE 


COEFFICIENTS 


V.  COMPUTATIONAL  ASPECTS 


A.  USING  THE  NUMERICAL  SOLUTION 

For  most  of  the  numerical  work,  the  state  probabilities  have  been 
obtained  using  the  fourth-order  RUNGE-KUTTA  method.  The  accuracy  of  the 
results  was  increased  by  substituting  the  available  analytical  results 
for  the  Regions  I,  II  and  III. 

The  disadvantage  is  that  this  solution  method  needs  a  lot  of  CPU-time. 
For  the  battle  in  Example  1,  with  the  data  given  in  Table  1,  the  calcula¬ 
tion  of  the  state  probabilities,  expected  force  levels  and  variances  led 
to  a  CPU-time  on  an  IBM-360  computer  of  almost  90  minutes  with  a  time 
step  size  of  0.05  minutes.  That  was  the  reason  why,  in  the  examples  the 
attrition  rate  coefficients  were  increased  by  one  magnitude,  which  brought 
the  computer  usage  down  to  a  CPU-time  of  around  12  minutes  with  the  same 
time  step  size  of  0.05  minutes.  The  sum  of  the  state  probabilities  at 
every  time  step  was  used  as  a  measure  of  accuracy.  Surprisingly,  its 
deviation  from  1.0  was  always  in  the  fourth  or  fifth  decimal,  which  proves 
the  robustness  of  the  RUNGE-KUTTA  method.  Also,  the  fact  that  for  the 
Regions  I,  II  and  III  analytical  solutions  were  substituted  at  each  time 
step  did  not  change  the  final  outcome  considerably.  It  was  found  that 
without  the  analytical  partial  solutions  the  sum  of  probabilities  tended 
to  be  slightly  higher,  changing  the  sum  of  probabilities  from,  for 
example,  0.99998  to  1.00003. 
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B.  THE  EQUAL  ATTRITION  RATE  COEFFICIENT  SOLUTION 

After  the  development  of  the  equal  attrition  rate  coefficient  solution 
(EARCS)  outlined  in  Ch.  III.C,  it  was  implemented  on  an  IBM-360  comput¬ 
er.  The  CPU-time  for  the  calculation  of  the  state  probabilities,  expected 
force  levels  and  variances  for  a  given  point  in  time  always  stayed  below 
20  seconds,  which  emphasizes  its  computational  advantage. 

Further  research  showed  two  shortcomings  of  the  EARCS,  which  are 
easily  overlooked.  The  coefficients  C(m,n)  vary  over  a  wide  range  start¬ 
ing  at  1.0  and  depending  on  the  initial  force  levels. 

For  example  for  a  battle  with 


mQ  =  20 


n0  =  40 


f  =  0.08 


%  =  15 

"bp  “  20 
tf  =  1.60995 


C(20 ,40)  =  1.0  but  C(16,21  )  =  4. 9297077276x1 0JJ. 

Also,  the  binomial  coefficients  in  equations  (3.23)  and  3.24)  vary  over  a 
wide  range  starting  at  1.0  to  1352078  for  the  above  example.  This  shows 
that  the  capacity  the  computer,  in  terms  of  the  number  of  significant 
digits,  is  able  to  carry  limits  the  implementation  of  this  solution  in 
the  present  form. 

Another  reason  why  this  solution  is  not  the  end  of  the  numerical 
problems  is  the  fact  that  in  the  given  form,  equations  (3.23)  and  (3.24) 
require  computation  of  an  alternate  sum  consisting  of  terms  whose  factors 
are  relatively  large  binomial  coefficients  and  small  numbers  between  zero 
and  one.  This  produces  truncation  errors,  which  yield  nonsensical  results 
like  negative  variances  and  sums  of  probabilities  greater  than  one. 
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VI.  CONCLUSIONS 


Throughout  this  thesis,  only  LANCHESTER-type  square-law  attrition 
with  fixed  initial  force  levels  and  fixed  attrition  rate  coefficients 
has  been  considered  for  a  deterministic  and  a  stochastic  version.  The 
stochastic  version  required  much  more  computational  effort.  So,  given 
the  need  for  an  analytical  model  as  opposed  to  the  use  of  simulation, 
there  is  not  much  to  gain  from  the  application  of  a  stochastic  model 
when  the  force  levels  are  large  and  the  forces  are  not  near  parity.  In 
these  cases  the  deterministic  version  essentially  produces  the  same 
results,  at  least  in  qualitative  terms.  For  smaller  force  levels  or 
forces  near  parity,  the  stochastic  version  may  be  helpful  to  get  more 
information  about  the  dynamics  of  combat. 

Also,  there  exist  one  case  where  the  deterministic  version  cannot 
differentiate  between  several  different  battles.  That  is  the  case  of 
equal  initial  forces,  equal  attrition  rate  coefficients,  but  different 
breakpoint  force  levels.  This  also  rectifies  the  further  development  of 
the  equal-attrition-rate-coefficient-solution  (EARCS). 

Using  an  analytical  model,  the  discussed  way  of  introducing  random¬ 
ness  into  the  model  seems  not  to  be  very  enlightening.  Therefore,  it  is 
suggested  that  another  way  to  include  random  effects  should  be  explored. 
The  author's  opinion  is  that  working  with  attrition  rate  coefficients 
which  are  random  variables  seems  more  promissing  to  gain  insight  into 
the  dynamics  of  combat.  Further  down  the  line  there  should  be  some 
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consideration  on  the  usage  of  combinations  of  the  possibilities  to 
include  random  effects  as  given  in  Chapter  II. B. 
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APPENDIX  A 


THIS  PROGRAM  COMPUTES  THE  STATE  PROBABILITIES  USING  THE  FOURTH 
ORDER  RUNGE-KUTTA  METHOD.  PLOTS  OF  THE  EXPECTED  FORCE  LEVELS  AND 


VARIANCES  AS  HELL  AS 
OBTAINED. 

OIHENSION  PI  133. 33J , 
DIMENSION  XMT  (500) . 

I  VMB  (5511  ,  VNBI55I1  , 


3-D-PLOTS  FOR  THE  STATE  PROBABILITIES  ARE 


P2I33.33),  T3DI20) 

XNT  (500)  .  OT  (500)  .  VM  (500)  .  VN  (5001  . 

VO  (55 1 )  ,  VARXI551).  VART(55I1,  TIM  (551 ) 


01  MENS  ION  EST  (3131  .0  (313)  .F  (2)  .SIZE  (2)  .  KX  (100)  ,  NT  ( 1001  ,NK  HI  I .  4 1 . 3) 


10 

1 

1  ,  PC  (41 ,41)  ,CST  (411  .CETERA  I 

11 

LOGICAL*!  10N(41.41I 

12 

DIMENSION  DETERM (313) 

13 

REALxB  TTL  (12)  /I2»’  V 

14 

REAL  Ml.  K2.  K3,  Kit 

IS 

CALL  ERRSET  (208,600.-1.  1) 

16 

F  (11-0. 

17 

FC2I-0. 

16 

LINES-0 

19 

SIZE  (1)  -6. 

20 

SIZE (2) -8. 

21 

OFLAG  •  0. 

22 

NKXT-100 

23 

READ  (5,21)  MBP,  MO.  NBP.  NO 

24 

READ  (5.221  AA.  BB 

25 

RE  AO  (5.22)  H.  FT  I M 

26 

READ (5,21)  N30 

27 

REAO (5.22)  EPSO.  EPSPTT 

28 

T30  (1)  •  9999. 

29 

IF  (N30  .LT.  11  GO  TO  25 

30 

REAO  (S, 22)  (T30  ( 1 ) ,  I-1.N3D) 

31 

25 

CONTINUE 

32 

C 

S3 

c 

N30-NUMBER  OF  30-PLOTS,  H-SIZE  OF  TIME  STEP 

34 

c 

FTIME-F1NALTIME.  EPSO.EPSPTT  ARE  ZERO  LEVELS 

35 

c 

T3D  (1) ■  FRACTION  OF  FT1ME  HHEN  TO  PLOT  3D 

36 

c 

37 

00  10  1-1. N3D 

36 

T30  (l)  -FT IM-T30  (I ) 

39 

10 

CONTINUE 

40 

21 

FORMAT  (1615) 

41 

22 

FORMAT  (8F10.5) 

42 

HR1TE  (6,28)  MBP,  MO.  NBP.  NO 

43 

28 

FORMAT!//  5X,  ’M  S  «  N  S'.  5X.  416  /) 

44 

HRITE  (6,29)  AA.  BB 

4S 

29 

FORMAT!/  SX,  'A  t  B',  9X.  2FI0.3  /) 

46 

HRITE  (6,30)  H.  FT1M 

47 

30 

FORMAT!/  SX.  'INITIAL  H  «  FINAL  TIME  (•  LOOPS 

46 

HRITE  (6,31)  N30 

49 

31 

FORMAT!//  SX.  '•  OF  3-D  PLOTS',  2X.  16  /I 

SO 

IF  (N30  . LT .  II  GO  TO  36 

•<  H)  \  SX.  2FI0.3/) 
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51 

WRITE  (6,241  EPSQ.EPSPTT 

52 

24  FORMAT  (/  SX.  ,EPSQ-\F10.5.5X. 'EPSPTT 

53 

WRITE  (6,  32)  (T3D  ( 1 )  ,  1-1, N30) 

54 

32  FORMAT  1/  5X.  * AT  TIME*.  2X,  10F10.3 

55 

36  CONTINUE 

56 

HOI  «  HO  ♦  1 

57 

NO)  -  NO  ♦  1 

58 

NS  -  N01  '  N8P 

59 

HS  •  HOI  -  HBP 

60 

HSL  «  HS  ♦  1 

61 

NSL  ■  NS  +  1 

62 

BMO  •  8B  m  HO 

63 

AND  -  AA  «  NO 

64 

OANB  -  ANO  /  BB 

65 

OBMA  •  BHO  /  AA 

66 

BAT  »  -  (BMO  ♦  ANO) 

67 

00  47  HI  «  l.HS 

66 

VM  (Ml)  ■  MBP  ♦  Ml  -  1 

69 

47  CONTINUE 

70 

00  48  N1  •  1 .NS 

71 

VN(Nl)  -  NBP  ♦  N1  -  1 

72 

48  CONTINUE 

73 

L  -  0 

74 

N3  -  1 

75 

TPLOT  «  T30  (N3) 

76 

TIME  •  0. 

77 

50  CONTINUE 

78 

IF  (TIHE  .GT.  FTlM)  GO  TO  210 

79 

52  CONTINUE 

60 

TIME  •  H  x  L 

81 

l  -  l  ♦  1 

82 

LB  -  0 

83 

TIM(L)  «  TIME 

84 

IF  (  L  .EQ.  21  GO  TO  53 

85 

IF  (TIME  .IT.  TPLOT)  GO  TO  60 

86 

C 

87 

C 

SET  NEXT  PLOT  TIHE 

88 

C 

89 

N3  •  N3  *  1 

90 

TPLOT  *  T30  (N3) 

91 

IF (N3.GT.N30!  TPUOT-9999. 

92 

53  CONTINUE 

93 

C 

94 

C 

LB  •  L  WHEN  TIME  TO  PLOT 

95 

C 

98 

LB  -  L 

97 

EBT  •  EXP  (BB  x  TIHE)  -  1. 

98 

EAT  -  EXP  (AA  «  TIME)  -  1. 

99 

AE  •  OANB  «  EBT 

100 

BE  •  OBMA  «  EAT 
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101 

BAE  •  EXP  (BAT  «  TIME) 

102 

fj  •  t. 

103 

JJ  -  1 

104 

C 

10S 

C 

NO  VALUE  OF  J  •  ZERO  AT  HO, NO 

106 

C 

107 

FH  ■  1. 

108 

KK  -  1 

109 

OFK-1. 

110 

OFJ-1. 

111 

C 

112 

C 

NO  VALUE  OF  K  *  ZERO  AT  HO. NO 

113 

C 

114 

60 

CONTINUE 

115 

IF  (L  . EO.  1)  GO  TO  500 

116 

00*0. 

117 

00  400  N!  •  I.N5 

118 

N  «  N01  -  N 1 

119 

NL  ■  NSL  -  NI 

120 

NL1  •  NL  ♦  1 

121 

ANJ  «  AA  m  N 

122 

C 

123 

C 

124 

00  300  HI  •  1 . MS 

12S 

H  •  HOI  -  HI 

126 

HL  •  MSL  -  Ml 

127 

ML1  «  ML  ♦  1 

128 

BMI  -  BB  «  M 

129 

ABA  «  ANJ  *  8H1 

130 

IF  IM  .EO.  MO)  GO  TO  502 

131 

IF  (M  .Ed.  M8P)  GO  TO  504 

132 

IF  (N  .EO.  NO)  GO  TO  506 

133 

IF  (N  .EO.  NBP)  GO  TO  507 

134 

C 

135 

C 

OEFAULTS  TO  ALL  INSIDE  POINTS 

136 

C 

137 

Kl  -  AN J«P2 (ML  1 .  NL)  ♦  BMI *P2 (ML ,  NL l )  -  ABA-P2 (ML.NL) 

136 

PT  -  ANJxO. 5*  (PMHLl.NLMP2tMLl.ND)  ♦  BMIxO.S  «  (P I  (ML , NL I )  ♦ 

139 

1  P2  (HL.NL1)  1 

140 

K2  »  PT  -  ABA  w  (P2IML.NL)  ♦  H«0.5«K1) 

141 

K3  ■  PT  -  ABA  «  IP2  (ML,  NL)  ♦  H«0.5«K2I 

142 

K4  -  ANJxPl  (HL1.NL)  ♦  BHI*P1  (ML, NL 1 )  -  ABA«  (P2  (ML .  NL <  ♦  H«K3) 

143 

PI  IHL.NLI  •  P2lML.NL)  ♦  (H/6.0)  «  (KM2.0«K2  ♦  2.0«K3  ♦  K4) 

144 

GO  TO  200 

145 

502 

IF  (N  .EO.  NO)  GO  TO  503 

146 

IF  (N  .EO.  NBP)  GO  TO  507 

147 

C 

148 

c 

H-MO.  N— NO. NBP 

149 

c 

150 

Kl  -  BH1xP2 (ML.NL11  -  RBA«P2 (ML.NL) 
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tSl  PT  •  611!  mO,  S  «  (Pl(ML.NLl)  ♦  P2IML.NL1II 

IS?  K2  -  PT  -  RBR  «  (P2lMt.NL)  ♦  HxkJxO.5) 

153  K3  ■  PT  -  RBR  x  (P21ML.NU  ♦  HxK2x0.5) 

154  K4  «  BN!  x  Pl(ML.NLl)  -  RBR  x  (P2IML.NU  ♦  HxK3> 

155  PI  (ML.  NL)  »  P2  (ML .  NL)  ♦  (H/6.)  x  (Kl»2.0xK2  ♦  2.0mK3  ♦  Km 

15S  if  (LB  .£Q.  U  GO  TO  601 

157  GO  TO  200 

158  C 

159  C  M«MO.  N»NO 

160  C 

161  503  CONTINUE 

162  Pl(ML.NL)  »  EXP  (-RBR  x  TIME) 

163  WRITE  (6, 102)  L.  N.  ML, NL .  Pl(HL.NL) 

161  GO  TO  200 

165  504  IF  (N  .EQ.  NBP)  GO  TO  505 

166  C 

167  C 

168  C  M-M6P.  N— «  NBP 

169  C 

170  K1  ■  RNJ  x  P2 (ML  1 , NL) 

171  K2  ■  RNJ  x  0.5  x  (PI(MLI.NL)  ♦  P2(MLI,NL)) 

172  K3  •  K2 

173  K4  «  RNJ  x  PI (ML1 .NL) 

174  PllML.NL)  -  P2  (ML , NL)  ♦  (H/6.0)  x  (KI  ♦  2.0xK2  ♦  2.0xK3  ♦  K4) 

175  GO  TO  200 

176  C 

177  C  M«HBP.  N*N8P 

178  C 

179  505  CONTINUE 

180  PI  (ML , NL)  •  0.0 

161  GO  TO  200 

182  C 

183  C  N-NO,  H— HBP.  MO 

184  C 

185  506  Kl  ■  RNJxP2(ML1,NL)  -  RBflxP2IML.NL) 

186  PT  -  RNJ  x  0.5  x  (PI  (ML  1 , NL)  ♦  P2iMLI.NL)) 

187  K2  •  PT  -  RBR  x  (P2(HL.NL)  ♦  Hx0.5xKI) 

188  K3  •  PT  -  RBR  x  (P2(ML,NL)  +  Hx0.5xK2) 

189  K4  •  RNJ  X  Pl(MLl.NL)  -  RBR  x  IP2IML.NU  ♦  HxK3) 

190  PllML.NL)  •  P2IML.NL)  ♦  (H/6.0)  x  (Kl  ♦  2.0xK2  ♦  2.0xK3  ♦  K4) 

191  IF  (LB  .EO.  L)  GO  TO  602 

192  GO  TO  200 

193  C 

194  C  N»NBP.  *  N*NBP.  M-MO 

195  C 

196  507  CONTINUE 

197  Kl  -  8MI  x  P2  (ML, NL 1 ) 

198  K2  •  BM1  x  0.5  x  (PI  (ML.  NL  1 )  ♦  P2  (ML ,  NL  1 ) ) 

199  K3  ■  K2 

200  K4  •  BHI  x  PI  (ML.NL1I 
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SOI 

SOS 

S03 

SOU 

505 

506 
307 

soe 

309 

310 
Sll 
SIS 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 


PI  (Ml. Ml)  ■  P2(ML,NLI  ♦  (H/6.0)  «  (K 1  ♦  2.0-K2  ♦  2.0-K3 
GO  TO  200 
C 

C  ALSO  M=MO.  N— NBP.NO 

c 

601  CONTINUE 

IF  (  (BE.LT.  1.1  .AND.  (KK.GT.9)  )  GO  TO  255 

C 

C  COMPUTE  PTT  IF  3E  >  0  4  KK  <  10 
C 

PTT-BAE 

DO  250  1ND* 1 , KK 
XI  NO* 1  NO 
PTT-PTT-8E/XIN0 
250  CONTINUE 
GO  TO  2S0 
2S5  PTT-O. 

260  CONTINUE 
KK  ■  KK  *  1 

IF  ( (ABS  (PI  (HL.NU  -PTT)  l  .LE.EPSPTT)  GOTO  199 
GO  TO  700 
C 

C  M-MBP.  N*NO 

c 

602  CONTINUE 

IF  1  (AE.LT.  1.)  .AND.  UJ.GT.9))  GO  TO  270 
C 

c 

C  COMPUTE  PTT  IF  AE  »  0  \  JJ  <  10 
C 

PTT-BAE 

00  265  IN0-1.JJ 
XI  NO- 1  NO 
PTT-PTT-AE/XINO 
265  CONTINUE 
GO  TO  275 
270  PTT-O. 

275  CONTINUE 
JJ  -  JJ  ♦  1 

IF  IIABSlPl  (ML.NLJ  -PTT))  .  LE.  EP5PTT)  GO  TO  199 

700  CONTINUE 
C 

C  TIME  REDUCED  1/2 
C 

L  •  L  -  I 

WRITE  (6,701)  TIME.  M,  N.  H 

701  FORMAT!//  5X.  'H  VALUE  IS  REDUCED  BT  HALF  AT  t  TIME  -  \ 

1  *  M  -  \  IS.  ’  N  -  *.  15,  ’  FROM  ’.  F8.3  //) 

H  ■  0.5  «  H 

WRITE (6, 102)  L.M.N.KK,  PTT.  PI  !ML . NL) . QFK.  QFJ.QO.OFLAG 


♦  K4) 


F8.3. 
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251  102  FORMAT (/  2X.  HIS,  7F12.5) 

252  QFLAG-O. 

253  GO  TO  52 

2511  199  PTT-P1  IML.NL) 

255  200  CONTINUE 

256  0Q-QQ*P1  IML.NL) 

257  300  CONTINUE 

258  400  CONTINUE 

259  OOJ-ABS  (00-1 .1 

260  IF (QQl . LE.EP5QI  GO  TO  401 

261  OFLAG*  1 . 

262  GO  TO  700 

263  401  CONTINUE 

264  C 

265  C  COHPUTE  VAfiX(T),  YARY  (T)  .  NBAR  IT  >  .  MBAR  (T)  .0  (T) 

266  C 

267  SMB  .  0. 

268  SN8  •  0. 

269  SO  -  0. 

270  SM2  -  0. 

271  SN2  -  0. 

272  DO  415  Ml  *  1,MS 

273  RM  «  VMIM11 

274  RM2  -  RM«RM 

275  00  410  N 1  «  l.NS 

276  RN  «  VN  INI I 

277  RN2  •  RNxRN 

278  PT  •  Pi  (M1.N1) 

279  P2IMI.N1I  »  PT 

260  SMB  -  SMB  ♦  RM«PT 

261  SNB  «  SNB  ■»  RN«PT 

282  SO  »  SO  ♦  PT 

283  SH2  •  SM2  ♦  RM2  «  PT 

264  SN2  •  SN2  ♦  RN2  «  PT 

285  410  CONTINUE 

286  415  CONTINUE 

287  VMS  III  »  SMB 

288  VN8IU  «  SNB 

289  VO  (LI  •  SQ 

290  VARX  (LI  ■  SM2  -  SMBxSMB 

291  VARY  (L)  *  SN2  -  SNB«.SNB 

292  C 

293  C  RETURN  TO  MAIN  LOOP  (50)  IF  NOT  TIME  TO  PRINT 

294  C 

295  IF  (N30  .LT.  II  GO  TO  430 

296  IF  (LB.  NE. LI  GO  TO  430 

297  IF  (LB. EQ. 21  GO  TO  430 

298  C 

299  C  ADJUST  X.Y.  VECTORS  FOR  PLOT 

300  C 
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30 1 

DO  15  K« 1 , MS 

302 

15 

EST  (K)  -K-l 

303 

00  16  K-I.NS 

304 

DETERM  IK)  »K-1 

305 

16 

CONTINUE 

306 

DO  115  1-1,41 

30? 

CST  (11-1 

308 

CETERM  (11-1 

309 

DO  115  J-1.41 

310 

PC(I,JI*0.0 

311 

115 

CONTINUE 

312 

MP-4I 

313 

NP*  41 

314 

flLP-15. 

315 

BETfl-30. 

316 

PHAX-P1  (1,1) 

317 

DO  116  1-1, NS 

318 

WRITE  (6. 199)  I.P1  (1  .1)  ,P 

319 

199 

EORMfiT  t *  *,I5. 10X.2F15.5) 

320 

116 

CONTINUE 

321 

00  18  1-1. MS 

322 

00  18  J-l.NS 

323 

IF  (PI  (1,  Jl  .GT.PMPIX)  PMRX-I 

324 

18 

CONTINUE 

325 

CONST-8. O/PMflX 

326 

00  19  1-1. MS 

327 

N-I-M8P 

328 

DO  19  J-l.NS 

329 

N-J-NBP 

330 

Pi  (1 ,  Jl  -CONST-PI  (t.J) 

331 

PC  (N. MI-PI  (I.JI 

332 

19 

CONTINUE 

333 

CALL  PLT301  (CST. MP, CETERM 

334 

1KXT, LINES) 

335 

C 

336 

430 

CONTINUE 

337 

C 

338 

C 

339 

C 

PLOT  HERE 

340 

C 

RETURN  TO  MAIN  LOOP 

341 

C 

342 

GO  TO  50 

343 

500 

CONTINUE 

344 

00  501  1  •  1 , MS 

345 

DO  501  J  *  l.NS 

346 

PI  (I.  J)  *  0.0 

347 

501 

CONTINUE 

348 

PI  (MS, NS)  •  1.0 

349 

GO  TO  401 

350 

C 

,  1  i 


,  BETfi.F.TTL.SIZE.WK. ION.KX.KT.N 


l 
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351  C  FINAL  TIME  REACHED 

352  C 

353  210  CONTINUE 

354  WRITE  16.2111  TIME.  FT 1M,  L 

355  211  FORMAT!//  5X,  ‘COMPUTED  TIME  \  F10.3.  5X.  ’INPUT  FINAL  TIME 

356  t  F10.3.  5X.  ’•  OF  LOOPS  TO  REACH  FINAL  TIME  \  15  //) 

357  00  215  I  »1 ,L 

358  WRITE  (6.213)  VMBlll.  VNBI1),  VQU),  VARXU).  VART(l),  TIM(I) 

35S  2l3  FORMAT I2X,  12F10.SI 

360  215  CONTINUE 

361  650  FORMAT  Cl  *1 

362  WRITE  (6.650) 

363  CALL  PLOTP  (TIM. VMB.L.O) 

36H  WRITE  (6,650) 

365  CALL  PLOTP(TIM.VARX.L.O) 

366  WRITE  (6.650) 

367  CALL  PLOTP  (TIM,  VNB.L.O) 

368  WRITE (6.650) 

369  CALL  PLOTP(TIM.VART.L.O) 

370  STOP 

371  DEBUG  SUBCHK 

372  END 
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APPENDIX  B 


1  c 

2  C  THIS  PROGRAM  CALCULATES  THE  STATE  PROBABILITIES  FOR  THE  EQUAL 


3 

C 

ATTRITION  RATE  COEFFICIENT  SOLUTION  lEPRCSl  . 

9 

C 

A 

3-0  PLOT  IS  PRODUCED  USING  THE  VERSATEC  PLOTTER, 

5 

C 

6 

IMPLICIT  REAL*8  (A-H.O-Z) 

7 

CALL  ERRSET  (208.2SS. 10. 11 

a 

DIMENSION  C  (50,501  .  M  (50)  ,  N  (50)  ,F  (50)  ,  PTMN  (50.50) 

g 

REAL*4  SIZE  (2)  ,FL  (2)  .  HK  (91 .91 . 3)  .  X  (91)  ,  T  (91)  .P  (9 

10 

DIMENSION  KX  (100)  ,  K t  (100) 

n 

LOGICAL*!  I ON  (4  1 , 9  1 ) 

12 

FiE AL*8  TTL  (121  / 1 2*>  *  V 

13 

RE AO  (5. 100) MO.MBP.NO. NBP 

19 

READ  (5. 101) A 

15 

READ  15. 102) TIME 

16 

READ  (5,  103)  IF  111  .  1.1,91) 

17 

c 

ie 

c 

F  (II -1/I-FACT0RIAL.  DONE  TO  SPEED  UP  THE  PROGRAM 

19 

c 

20 

WRITE  (6,802) MO.MBP.NO. NBP.  A.  TIME 

21 

602 

FORMAT  ( 1  ’ , 4 1 5. 2F 10. 5) 

22 

MO-M0-M8P 

23 

NO-NO-NBP 

29 

M01-MO-M 

25 

NO  1  -NO*  1 

26 

MBP1-MBP-M 

27 

NBP1=NBP»1 

20 

M0P2-MBP»2 

29 

N8P2»NBP*2 

30 

RMO-MO 

31 

RNQ-NO 

32 

00  10  1*1,41 

33 

DO  10  J.1,41 

39 

C<!,J)*0.0 

35 

PTMN  (I.  Jl  *0.0 

36 

Pll.Jl-O.O 

37 

10 

CONTINUE 

38 

C(MOl.NOl)  =1.0 

39 

00  20  1*2. MO 

90 

J*MBP»I 

91 

C  (J.NOl)  «  RNO«*  (MO-  I-  1) 

92 

20 

CONTINUE 

93 

DO  21  1-2.N0 

99 

J-NBP-U 

95 

C  (M01 ,  Jt  •  RHO*«  (NO-  I  ♦  1 ) 

96 

21 

CONTINUE 

97 

DO  22  1*2, MO 

98 

MC-M01-I 

99 

MM*M01 -  I ♦ 1 

50 

MPLUS-MM*] 

81 


51 

DO  22  J*2.ND 

52 

NC*N01-J 

53 

NN*NQ  1  -  _)♦  1 

54 

NPUiS^NN+l 

55 

C(MH.NN)«NC«  C (HPLU5.NN) *MC«  C(MM.NPLUS) 

56 

22 

CONTINUE 

57 

C 

58 

C 

END 

OF  COEFFICIENT  CALCULATION 

59 

C 

60 

FI—  AxTIME 

61 

DO  25  1  =<MBP2,  MO  1 

62 

00  25  J=NBP2, NQ1 

63 

64 

K«MO-'NQ-IJ  +2 

65 

F2* (1 .0-DEXP  (Fill  —K 

66 

KK  =  K*  1 

67 

IJ*IJ~2 

68 

F3-0EXP  m«i  J1 

69 

PTMN  (1.  J1  *F2«F3xC  11.  J)  *F  IKK) 

70 

25 

CONTINUE 

71 

100 

FORMAT  (415) 

72 

101 

FORMAT  (F10.S) 

73 

102 

FORMAT  (F10.5) 

74 

103 

FORMAT  (017.  11) 

75 

c 

76 

c 

FOR 

nbp<n<no 

77 

c 

76 

DO  30  11-l.ND 

79 

NF0RCE*N8P*I 1 

60 

J*MO«NO-MBP-1-NFORCE 

81 

fact»mbp-m.o+nforce 

82 

SUM1  *  (1 .0-DEXP  (FI -FACT))  /  1  A«FACT) 

83 

5UMI NT  *SUM 1 

84 

00  301  K*1.J 

85 

FACT-FACT+l  .0 

86 

ADDFAC*  (-1.0)  «»K 

87 

FRACTN*  1 1 .  Q-DEXP (Fl-FOCT) ) )  t  A-FACT) 

88 

COMBT-I .0 

89 

DO  302  K J= 1 ,K 

90 

RKJ*KJ 

91 

COMBI*  (J-fiKJ-M  .  0)  /R HJ 

92 

C0MBT*C0H8T  «C0M8 1 

93 

302 

CONTINUE 

94 

SUM1 NT* SUM  I  NT ♦AODFAC»FRACTN-COMBT 

95 

301 

CONTINUE 

96 

NM8P* J+l 

97 

N80UN0*NF0RCE  *  1 

98 

PFAC*A«NFORCEkF 1NMBPI  -C (MBP2. NBOUND) 

99 

PTMN 1M8P1 . NC0UND1  =PFAC*SUM1NT 

100 

30 

continue 
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101 

C 

102 

C 

FOR  MBP<M<MO 

103 

C 

10>1 

00  40  11=1. MO 

105 

MF0RCE-M8P-1 1 

106 

J-MO-NO-NBP-l-MFORCE 

107 

FACT-NBPM  .O-HFORCE 

108 

SUM1*  (l.O-OEXP  (FI -FflCTn  /(  ft-FRCT) 

109 

SUH1NT=SUM1 

110 

DO  401  K=1.J 

111 

FRCT«FflCT ♦1.0 

112 

ADOF RC«  (-1.0) ««K 

113 

FRACTN-  (l.O-OEXP  (F  l-FACT)  )  /  l  A-FACT) 

114 

COMBT  =  1.0 

115 

00  402  K J* I , K 

116 

RK J*KJ 

117 

COMBI*  1 J-RK  J*  1.01  /RK  J 

118 

COMB  T  =  C0M8T><  COMB  I 

119 

402  CONTINUE 

120 

SUM 1  NT  =  SUM  1  NT ♦ROOF  AC-FRACTN-COMBT 

121 

401  CONTINUE 

122 

MNBP-J-1 

123 

M80UN0*MF  QRCE* 1 

124 

PFAC-A.MFORCE-F  (MNBP) «C  IMBOUHD , NBP2) 

125 

PTMN  (MBOUNO. NBPl 1 -PFAC-SUMINT 

126 

40  CONTINUE 

127 

00  45  I=MBPl.MOl 

128 

00  45  J=N8P 1 , NO  1 

129 

K-l-1 

130 

L-J-l 

131 

WRITE  (6.801IK.L.PTMNU.  J) 

132 

801  FORMAT  ( "  P(T.\  15. •.MS.’)  =  ’.017 

133 

45  CONTIN'  E 

134 

C 

135 

c 

&ATA  ADJUSTMENT  FOR  PLOT 

136 

c 

137 

NROM-41 

138 

NCOL-41 

139 

NKXT  *  1 00 

140 

UNES-0 

141 

ALPHA. 15. 

142 

BETA-30. 

143 

FL  (11-0.0 

144 

Fl  (21-0.0 

145 

SIZE  (1)  -6.0 

146 

SIZE  (21  -8.0 

147 

c 

148 

c 

SCALING 

149 

c 

150 

PMAX-PTMN  (1 . 1) 
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151  DO  50  1-1.41 

152  XU)- 1 

153  00  50  J-1,41 

154  TIJl-J 

155  IF  IPTHN  11.  J)  .GT.PMAX1  PMAX-PTMNlI,  J) 

156  SO  CONTINUE 

157  CONST-8. O/PMAX 

156  DO  51  1-1.41 

159  DO  51  J-1,41 

160  PII.JI  -  SNGL  (CONST -PTMN  11.  J) ) 

161  51  CONTINUE 

162  C 

163  C 

164  C  IF  OTHER  FORCE  LEVELS  CHANGE  PLOT  ARGUMENTS 

165  C 

166  CALL  PLT301  (X.NROH.T.NCOL, P. ALPHA. BETA. FL.TTL. SIZE. WK. 

167  «  lDN.KX.KT.NKXT. LINES) 

168  C 

169  C  EXPECTED  VALUES  AND  VARIANCES 

170  C 

1 7 1  EM-0.0 

172  EMM-0.0 

173  VARM-0.0 

174  SPR08-0. 0 

175  00  60  1I-MBP1.M01 

176  1-11-1 

177  SPT-0.0 

178  DO  601  JJ-NBP1.N01 

179  SPT-SPT  ♦  PTMNU1.JJI 

180  601  CONTINUE 

101  EM-EM* 1 wSPT 

182  EMM-EMM*  I-I-SPT 

183  SPROB-SPRO0*SPT 

184  60  CONTINUE 

185  VARM-EMM-EM-EM 

186  EN-0.0 

187  ENN-0.0 

188  VARN-0.0 

189  DO  70  11-N8P1.N01 

190  1-11-1 

191  SPT-0.0 

192  DO  701  JJ-MBPI .MO! 

193  SPT-SPT*  PTMNUJ.II) 

194  701  CONTINUE 

195  EN-EN* 1 -5PT 

196  ENN-ENN*  I-l-SPT 

197  70  CONTINUE 

198  VARN-ENN-EN-EN 

199  WRITE  16.805)  EM. EN. VARM. VARN 

200  805  FORMAi  <  *  '.’EXPECTED  VRLUE5  M  .  N\  SOX.’  VARIANCE  M.N 
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